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Abstract 

We construct the chiral effective Lagrangian for two lattice theories: one with Wilson fermions 
and the other with Wilson sea fermions and Ginsparg- Wilson valence fermions. For each of these 
theories we construct the Symanzik action through 0(a 2 ). The chiral Lagrangian is then derived, 
including terms of 0(a 2 ), which have not been calculated before. We find that there are only few 
new terms at this order. Corrections to existing coefficients in the continuum chiral Lagrangian are 
proportional to a 2 , and appear in the Lagrangian at 0(a 2 p 2 ) or higher. Similarly, 0(4) symmetry 
breaking terms enter the Symanzik action at 0(a 2 ), but contribute to the chiral Lagrangian at 
0(a 2 p ) or higher. We calculate the light meson masses in chiral perturbation theory for both 
lattice theories. At next-to- leading order, we find that there are no 0(a 2 ) corrections to the 
valence-valence meson mass in the mixed theory due to the enhanced chiral symmetry of the 
valence sector. 
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I. INTRODUCTION 

Chiral perturbation theory (x?T) plays an important role in the analysis of current 
lattice QCD data. Simulations with the quark masses as light as realized in nature are not 
feasible on present day computers. Instead one simulates with heavier quark masses and 
performs a chiral extrapolation to the physical quark masses using the analytic predictions 
of %PT. To perform the chiral extrapolation one must first take the continuum limit of 
the lattice data, since ;^PT describes continuum QCD and is not valid for non-zero lattice 
spacing. However, it is common practice not to perform the continuum extrapolation and 
nevertheless fit the lattice data to continuum ;^PT, assuming that the lattice artifacts are 
small. 



A strategy to reduce this systematic _uncertainty was proposed in Refs. |3j 
different approach was taken in Refs 



a (a 



• 00] in the strong coupling limit). There it was shown 
how the discretization effects stemming from a non-zero lattice spacing can be included 
in ;\RT. The basic idea is that lattice QCD is, close to the continuum limit, described 
w Symanzik's effective theory, which is QCD with additional higher dimensional terms 
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IkI ll^ . The derivation of x?T from QCD can then be extended to this effective 
theory with additional symmetry breaking parameters. The result is a chiral expansion in 
which the leading dep endence on the lattice spacing is explicit. This idea was numerically 



examined in Ref. 
the results of Ref. 



for a theory with two dynamical sea quarks on a coarse lattice using 
. The characteristic chiral log behavior in the pseudo scalar meson 
mass and decay constant was observed. 

A similar approach was taken in Ref. 16] for analyzing lattice theories with two types 
of lattice fermions - Wilson fermions for the sea quarks and Ginsparg- Wilson fermions 



for the valence quarks. The latter can be implemented using domain wall 



overlap 
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23, 



24j . perfect action 



25l |26| , and chirally improved fermions |2j 



19|, 



There are several advantages in using different lattice fermions in numerical simulations. 
Since massless Ginsparg- Wilson fermions exhibit an exact chiral symmetry even at non-zero 
lattice spacing , it is possible to simulate such valence fermions with masses much smaller 
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31]. 



This allows a wider 



than the valence quark masses accessible using Wilson fermions 
numerical sampling of points in the chiral regime of QCD. In addition, the valence sector 
exhibits all the benefits stemming from the Ginsparg- Wilson relation 32] such as the absence 



2 



of additive mass renormalization, of operator mixing amongst different chiral multiplets and 



!3U2J, I25U26 
;. H and |lj 



of lattice artifacts linear in the lattice spacing a 

In this paper we extend the results of both Refs. |7[ and |16( by calculating the chiral La- 
grangian including the 0(a 2 ) lattice effects. There are various reasons for doing this. First, 
the lattice spacings in current unquenched simulations are not very small, so that neglect- 
ing the 0(a 2 ) contributions might not be justified. Second, the use of non-perturbatively 
improved Wilson fermions in lattice simulations is becoming more common. The leading 
corrections for these fermions are of 0(a 2 ) and hence need to be computed in order to know 
how the continuum limit is approached. 

At 0(a 2 ) many operators enter the Symanzik action and need to be taken into account 
for constructing the chiral Lagrangian. 4-fermion operators appear for the first time and 
operators that explicitly break Euclidean rotational symmetry are encountered. Neverthe- 
less, the number of new operators in the chiral Lagrangian is rather small (3 for the Wilson 
action and 4 for the mixed fermion theory). This is important in practical applications, since 
every new operator comes with an undetermined low-energy constant. These constants en- 
ter the analytic expressions for physical observables and too many free parameters limit the 
predictability of the chiral extrapolations. 

The paper is organized as follows: %PT for the Wilson action is discussed in section [TT1 
including the partially quenched case in subsection III El The mixed theory with Wilson sea 
and Ginsparg- Wilson valence quarks is treated in section IIHI In section US we discuss the 
chiral power counting and compute the pseudo scalar meson mass including the new 0(a 2 ) 
contributions for both cases. We end with some general comments in section 



II. WILSON ACTION 



In this section we formulate the chiral effective theory for the Wilson lattice action. First 
the Wilson action and its symmetries are briefly reviewed, then the local Symanzik action 
through 0(a 2 ) is presented. Based on the symmetry properties of the Symanzik action we 
construct the chiral effective theory. Finally, we consider the extension to the partially 
quenched case. 
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A. Lattice action 



We consider an infinite hypercubic lattice with lattice spacing a. The quark and anti 



quark fields are represented by ip and t/j, respectively. Wilson's fermion action 
by 



is given 



S w = a A ^2?p(D w + m )ip(x) , (1) 

X 

Dw = l{l»(K + V»)-arV;Vj} , 

where mo denotes the Nf x Nf bare quark mass matrix and r the Wilson parameter. V* , V M 
are the usual covariant, nearest neighbor backward and forward difference operators. 

The Wilson action in eq. possesses several discrete symmetries - charge conjugation, 
parity - as well as an SU(N C ) color gauge symmetry. The introduction of a discrete space- 
time lattice reduces the rotation symmetry group 0(4) to the discrete hypercubic group. 

Next, we consider the group of chiral flavor transformations, 

G = SU(N f ) L <g> SU(N f ) R . (2) 

Introducing the usual projection operators P± = |(l±7s) the left- and right-handed fermion 
fields are defined by 

lpL,R = P^, ^ L>R =^P±. (3) 

Under a transformation L ® R G G these chiral components transform according to 

y LipL, i> L ^i) L L\ (4) 
tp R -> Rip R , ip R -> ip R R j . 

Only if the mass and the Wilson terms are absent is G a symmetry group of the theory. If 
all the quark masses are non-zero but equal the vector subgroup with L = R is a symmetry 
of the action Sw 

To complete the definition of the lattice theory one should also define a gauge action 
Sym- However, the precise choice of the gauge action is irrelevant for the purpose of our 
analysis, so we leave it unspecified. 
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B. Symanzik action 



The Symanzik action for the Wilson lattice action, up to and including 0(a 2 ), has been 
calculated first in Ref. [13]. The analysis to 0(a) has been later elaborated on in Ref. jl^ . 
We restate these results in a slightly different form. We list all operators in the action 
through 0(a 2 ) that are allowed by the symmetries, organized in powers of a (again, we only 
focus on the fermion action). We use the notation 



where O- are operators of dimension n constructed from the quark and gauge fields and 
their derivatives. The constants cf* are unknown coefficients. 

Some allowed operators in the Symanzik action are obtained by multiplying lower dimen- 
sional operators with powers of the quark mass m . Such operators (except for the mass 
term itself) renormalize lower dimensional operators. For example, a tr , (mo)^^ contributes 
to the wave function renormalization of the quark fields. Such operators can be effectively 
accounted for by an appropriate modification of the coefficients of the lower dimensional 
operators, and we do not list them explicitly. 

We find the following list of operators (we use the same notation as in Ref. Q|): 



Ss = a l S-i + S + aSi + a 2 S 2 . . . , 



(5) 




S. 



o • 



-l 



Of = tpip . 

of = ^V, 




(6) 
(7) 



By construction, Sq is the continuum QCD action. 



5- 



i 






5*2, bilinears: 



Of = v^V, 

Of = ^{D^D^U) + $D^D^ 



Of = ^D^D^, (9) 
Of = ^D^D^D^. 
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/~ T ±a i \ 2 

= {ipt^ipy, 


(6) 




(6) 

^13 




(6) 




0(6) 


= (^t°<vV) 2 , 



(10) 



where t a are the SU (N„ ) generators. This list of 4-quark operators is slightly different 
from the one in Ref. 13|- Sheikholeslami and Wohlert's list contains operators with flavor 
group generators. However, both lists are equivalent and are related by Fierz identities (see 
appendix lA"|) . Our choice of operators is guided by the fact that for the study of chiral 
transformation properties it is more convenient to consider 4-quark operators with trivial 
flavor structure. 

Equations of motion can be used to reduce the number of operators in the action. This 
also involves changing the definitions of the continuum fields that are matched to the lattice 
ones QQ. The way an operate, in the S ym an zik action affeets the ehna, L„n 
is essentially determined by its transformation properties under chiral rotations. By using 
equations of motion one can replace an operator only with operators in the same repre- 
sentation of the flavor symmetry group. Hence, the form of the chiral Lagrangian remains 
unchanged whether equations of motion are used or not. In the context of on-shell improve- 



HQ. 



ment, equations of motion have been used to reduce the number of operators at 0(a) 

(5) 

Only the Pauli term 2 is left at this order, and can be subsequently canceled by adding 
the clover term to the lattice action. For the formulation of the chiral effective theory we 
choose not to take this step. 

We distinguish two types of operators in the Symanzik action: those that break chiral 
symmetry and those that do not. Among the first type, the largest symmetry breaking 
effect comes from the term a _1 0[ 3 \ This term can be absorbed into by introducing the 
renormalized quark mass 

m = m q Z m , m q = m -m c , (11) 

where the critical quark mass m c ~ 1/a can be defined as the value of m where the pion 
becomes massless. Similarly, the effect of some higher dimensional operators can be absorbed 
in the definition of the renormalized quark mass by including a factor 1 + biam q + b 2 (am q ) 2 
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on the right hand side of eq. (fTTj) 



14. 



Q 

At 0(a) all operators break chiral symmetry. At 0(a 2 ) there are six symmetry breaking 
operators: Of\ Of\ Of \ Of$, 0$ and Of}. Fermionic operators that do not break the 
chiral symmetry first appear at 0(a 2 ). Purely gluonic operators (which we have not listed 
above) also belong to the second type of operators as they are trivially invariant under chiral 
transformations. They too enter at 0(ct 2 ). 

The operator Of deserves special attention. While respecting the chiral symmetries it 
is not invariant under 0(4) rotations. This means that it does affect the structure of the 
chiral Lagrangian by inducing 0(4) symmetry breaking terms in it. The analysis leading to 
the Symanzik action reveals that such terms must be at least of 0(a 2 ). 



C. Spurion analysis 

At 0(a°), the Symanzik action is QCD-like. For small a and m we assume the lattice 
theory to exhibit the same spontaneous symmetry breaking pattern SU(Nf)i^)SU(Nf)R — > 
SU (Nf) v as continuum QCD. Consequently, the low-energy physics is dominated by Nambu- 
Goldstone bosons which acquire small masses due to the soft explicit symmetry breaking 
by the small quark masses and discretization effects. The low-energy chiral effective field 
theory is written in terms of these light bosons. 

To construct the chiral Lagrangian we follow the standard procedure of spurion analysis. 
We write a term in the Symanzik Lagrangian as CqO where O contains the fields and their 
derivatives, and Cq is the remaining constant factor. For symmetry breaking terms, O 
changes under a chiral transformation of the fermionic fields, O — > O'. We then promote Co 
to the status of a spurion C, with the transformation C — > C such that CO = CO'. The 
chiral effective theory is constructed from the Nambu-Goldstone fields and the spurions with 
the requirement that the action is invariant under chiral transformations if the spurions are 
transformed as well. Once the terms in the chiral Lagrangian are obtained, each spurion 
is set to its original constant value C = Cq. This procedure guarantees that the chiral 
effective theory explicitly breaks chiral symmetry in the same manner as the underlying 
theory defined by the Symanzik action, and reproduces the same Ward identities. 

It might appear that one needs many spurion fields to accommodate all the symmetry 
breaking operators in the Symanzik action. However, this is not the case. Two spurions 
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which transform in the same way will lead to the same terms in the chiral Lagrangian, 
and so it is enough to consider only one of them. This is discussed in appendix [BJ Since 
we organize the chiral perturbation theory as an expansion in m and a, we do distinguish 
between spurions that transform the same way but have different m or a dependence. 

In the following we list the representative spurions. Shown are the transformation rules 
for the different spurions under chiral transformations, and the constant values to which the 
spurions are assigned in the end. 

0(a°) : M -> LMR\ M t -> RM^L\ (12) 

M = = m = diag(m 1; . . . ,m Nf ). 

This makes the mass term ip L M%l) R + iJ) R M'ip£ invariant under the chiral transformations of 
eq. (@J). 

0(a) : A -> LAR\ A ] 

A = Al = al, 

where / is the flavor identity matrix. 

0(a 2 ) : B = B 1 (g)B 2 ^ LB^ ® LB 2 R\ 5 f i 
C = C X ®C 2 ^ RCJ) ® LC 2 R\ C f e 
B = B\ = C = Cl = a 2 1 ® I. 

These spurions are introduced to make the symmetry breaking 4-quark operators invariant 
and therefore carry four flavor indices (see Ref. [37I and references therein). Consider, for 
example, the operator 

(^)WO = ($M($M + $r^l)$M + &rtPl)& r iPl) + $M$ R i>L). (is) 

The first term on the right hand side can be made invariant with the spurion B as can be 
seen from 

Bip L ip R ip L ip R = B ijk t(ip L )i(ip R ) j (iJ; L ) k (iJ; R ) l = tp L Bxip R tp L B 2 tfj R . (16) 

Similarly, all the other symmetry breaking 4-quark operators can be made invariant using 
the spurions B, C and their hermitian conjugates. At 0(a 2 ) one also needs to consider 
the squares of the 0(a) spurions. Note, however, that A 2 , A^A, (A^) 2 and AA* transform 
exactly like B, C, B^ and C>, respectively, and therefore need not be treated separately. 
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RA*L\ (13) 



: B\ ®B\^ RB\L ] <g> RB\l\ (14) 
C\ ® C\ -> LC\R) ® RC\L\ 



D. Chiral Lagrangian 



The chiral Lagrangian is expanded in powers of p 2 , m and a. Generalizing the stan- 
dard chiral power counting, the leading order Lagrangian contains the terms of Q(p 2 ,m,a), 
while the terms of 0(p 4 , p 2 m, p 2 a, m 2 , ma, a 2 ) are of next-to-leading order. In terms of the 
dimensionless expansion parameters mj A x and aA x , where A x « lGeV is the typical chiral 
symmetry breaking scale, this power counting assumes that the size of the chiral symmetry 
breaking due to the mass and the discretization effects are of comparable size. 1 

For the Wilson action, all next-to-leading order terms have already been computed in 
Ref. Q, except for the 0(a 2 ) terms. We are now in the position to calculate these contri- 
butions, which are the ones associated with the spurions B and C . We find the following 
three new terms (and their hermitian conjugates) 

(5 1 Et>( J B 2 Et>^a 2 (Et) 2 , (17) 
(B^B^) ^a 2 <EtEt), (18) 
<C 1 E)<C 2 Et)^a 2 <E><Et). (19) 

Here £ = exp(2zll//), with II being the matrix of Nambu-Goldstone fields. £ transforms 
under the chiral transformations in eq. (jlj) as £ — > LYiB) . The angled brackets are traces 
over flavor indices, and the arrows indicate assigning B = B , C = C Q , according to eq. (JT3J). 

So far we only considered the operators in the Symanzik action that explicitly break 
chiral symmetry. Operators that do not break chiral symmetry also contribute at 0(a 2 ). 
These operators do not add any new terms to the chiral Lagrangian, but simply modify the 
coefficients in front of already existing operators. At leading order, for example, the kinetic 
term is «4- (c^Ec^E^. There are corrections to f 2 due to the symmetry conserving terms 
in the Symanzik action: f 2 ^f 2 + a 2 K (K is another unknown low-energy constant.) This 
leads to the correction a 2 K (c^EdEj,) for the kinetic term. Thus given a term of 0(p 2 ) 
there is another term of 0(a 2 p 2 ). In general, we can rewrite the coefficient of any allowed 
operator in the chiral Lagrangian to obtain a new allowed operator which is 0(a 2 ) higher. 
These terms are beyond next-to-leading order and are not included in the present work. 

As already mentioned, the operator O4 breaks the 0(4) symmetry in the Symanzik 
action. However, in order to break the 0(4) symmetry, while still preserving the discrete 

1 A more detailed discussion of the power counting scheme is given in section Hvl 



9 



hypercubic symmetry, an operator must carry at least four space-time indices. In the chiral 
Lagrangian, these are provided by the partial derivative d^, hence the operator is at least 
of 0(p A ). Adding the fact that it is also an 0(a 2 ) effect, we see that the leading 0(4) 
symmetry breaking terms in the chiral Lagrangian are of 0(p A a 2 ) (an example is the operator 
a 2 (d^d^d^d^T,^). Hence, up to the order considered here, 0(4) breaking terms can be 
excluded from the analysis. 

Finally we can write down the terms of 0(a 2 ) which enter the next-to- leading order chiral 
Lagrangian. In terms of the two parameters 

m = 2B m = 2B diag(mi, . . . , rriN f ), a = 2Woa, (20) 

which have been introduced in Ref . 2 these terms are 

£ [a 2 ] = -a 2 W's{Y) + E) 2 - a 2 W 7 (^ - E) 2 - a 2 W' & {^ + EE). (21) 

The coefficients W- are new unknown low-energy constants. Putting it all together, also 
quoting the terms in the Lagrangian of 0(a) from Ref. j^, 3 we find 

£ x =y (^Ec^E 1 ") - ^ {rfiT) + Em) - <E f + E) (22) 

- L x <^E^Et) 2 - L 2 {dpEdvTt) (0 M £$,£ + ) - L 3 ((^Ec^) 2 ) 
+ U (d^d^) (mE f + Em) + aW A (d^d^) (£ f + E) 

+ L 5 ((9 M E(9 M E t (mE t + Em)) + aW 5 (d^d^(tf + E)) 

- Lq (mE f + Em) 2 - aW 6 (mS f + Em) <E f + E) 

- L 7 <mE t - Em) 2 - aW 7 (mE f - Em) <E f - E) 

- L$ (mTrmTr + EmEm) — dWg (wiY^Y) + EEm) 
+ L [a 2 ] + higher order terms. 

Here, the parameters Lj are the usual Gasser-Leutwyler coefficients of continuum xPT. 



2 Unlike in Ref. Q|, here we define d without the factor of c$w- The coefficient csw is not kept explicit 

(5) (5) (5) 

as we do not use equations of motion, and Si contains ac\ CL besides the Pauli term acs\y0 2 ■ 

3 There are some typos in the Lagrangian [eq. (2.10)] in Ref. pfl. These are corrected in eq. (J52J. 
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E. Partially quenched QCD 



Part ially quenched QCD is formally represented by an action with sea, valence and ghost 
quarks [38[. We collect the quark fields in \l> = (ips, ipy), where ips describes the sea quarks, 
and ipy contains both the anticommuting valence quarks and commuting ghost fields. The 
same is done for the anti quark fields. The mass matrix is given by m = diag(m,s, m.y), 
with ms being the Nf x Nf mass matrix for the sea quarks and m v = diag(my, my) is the 
2Ny x 2Ny mass matrix for the valence quarks and valence ghosts. 

We consider partially quenched lattice QCD with Wilson's fermion action eq. (0) for 
all three types of fields. The discrete symmetries and the color gauge symmetry is as in 
the unquenched case. The group of chiral flavor transformations, however, is different. If 
all the masses and the Wilson parameter r are set to zero, the action is invariant under 
transformations in the graded group 4 

G PQ = SU(N f + N V \N V ) L ® SU(Nf + N V \N V ) R . (23) 

Based on the symmetries of the lattice theory the Symanzik action for partially quenched 
lattice QCD is obtained as before. The result is easily quoted: One can simply replace ip 
and i/j in the Symanzik action for the unquenched theory with the extended fields \I/ and 
\I/ because the only two- and four-quark operators that are invariant under the extended, 
graded flavor group are still and its square. 

The leading term in the Symanzik action is partially quenched QCD, for which the 
construction of the chiral Lagrangian (first introduced in Ref. Q|) is essentially the same 
as for the unquenched case j39|. This remains true when higher dimensional operators in 
the Symanzik action are included, and the analysis of the previous subsection is readily 
extended to the partially quenched case. In particular, the form of the chiral Lagrangian 
for partially quenched lattice QCD with Wilson fermions is exactly the same as in eq. (|22j) . 
The difference is in the definition of the angled brackets, which now denote super-traces, 
and the interpretation of E and m. These need to be appropriately redefined to reflect the 
larger flavor content of partially quenched ;\TT. 



See Ref. for a more honest discussion of the symmetry group of partially quenched QCD. 
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III. MIXED ACTION 



In this section we consider a lattice theory with Wilson sea quarks and Ginsparg- Wilson 
valence quarks. As before we first construct Symanzik's effective action through 0(a 2 ). We 
then derive the chiral Lagrangian for this theory. 



A. Lattice action 



The use of different lattice fermions for sea and valence quarks is a generalization of 
partially quenched lattice QCD. Theoretically it too is formulated by an action with sea and 
valence quarks and valence ghosts. However, in addition to allowing different quark masses 
(mg 7^ my), the Dirac operator in the sea sector is chosen to be different from the one for 
the valence quarks and ghosts. For this reason we will refer to this type of lattice theories 
as "mixed action" theories. 

The mixed action theory with Wilson sea quarks and Ginsparg- Wilson valence quarks is 
defined in Ref. [16]. We refer the reader to this reference for details and notation. Here we 
just quote that the flavor symmetry group of the mixed lattice action is 

G M = G Sea ® G Val , (24) 
G Sea = SU{N f ) L <g> SU(Nf) R , 

G Val = SU(N V \N V ) L ® SU(N V \N V ) R . 

The quark mass term in the mixed action breaks both Gs e a and Gy a i- However, in the 
massless case Gy a i becomes an exact symmetry 2^] while Gs ea is still broken by the Wilson 
term. Because of the different Dirac operators there is no symmetry transformation that 
mixes the valence and sea sectors, in contrast to the partially quenched case [cf. eq. (J23J)]. 



B. Symanzik action 

The Symanzik action for the mixed theory can be derived using the results of the previous 
section. It is convenient to separately discuss three types of terms - those that contain only 
sea quark fields, those that contain only valence fields, and those that contain both. 

For the first type of terms the analogy with the previous section is evident: the relevant 
symmetry group is Gs e a — G, and the explicit symmetry breaking structure is the same. 
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Thus, all bilinear operators (if)) and 4-quark operators o[ (ip,ip), listed in subsection 
III Bl appear in Symanzik's action, once if) is replaced by ips- 5 

The construction of the purely valence terms is also analogous to the one for the Wilson 
action in subsection III Bl However, there are stricter symmetry constraints for Ginsparg- 
Wilson quarks and ghosts because the Ginsparg- Wilson action possesses an exact chiral 
symmetry when the quark mass is set to zero. All operators without any insertions of 
the quark mass must therefore be chirally invariant. In particular, all dimension-3 and 
dimension-5 operators are forbidden. Several dimension-6 operators at 0(a 2 ) are also ex- 
cluded. Only the bilinears of eq. (jHJ) and the 4-quark operators Of\ipv, Vv), i — 8, 12, 
and 13, of eq. (fTU|) . are Gy a i invariant and are therefore allowed. 

For terms of the third type, note that the symmetry group Gm forbids bilinears that mix 
valence and sea quarks. Thus, the only terms containing both sea and valence fields are 
4-quark operators which are products of two bilinears - one from each sector. Again, only 
the four terms Of\if>s, Vv), i — 7, 8, 12, and 13, are allowed. All the others break the 
chiral symmetry in the valence sector when my = 0. 

From these considerations it follows that the Symanzik action for the mixed lattice action 
up to and including 0(a 2 ) contains the following terms: 

S-i : Of^s), (25) 

So: Of\iP s ), Cf(Vv), i = l,2. (26) 

Si: 0l 5 \if> s ), i = l,2. (27) 

S 2 , bilinears : Of\ifj s ), Of\ifj v ), i = 1 - 4. (28) 

Sa, 4-quark operators : Of (tps, ips), i = 5 — 14, (29) 

OfO/v^v), Of\ip s ,iPv), i = 7,8,12,13. 

As before we have not listed operators that are obtained by multiplying lower dimensional 
operators with powers of the bare mass. 



5 We make the dependence of bilinear operators on the fields explicit by writing 0(ip). All the 4-quark oper- 
ators that we consider have the structure 0(ipi, V^) = ^ 1 il' 7 V ; iV ; 2^' / ' ( /'2- Here ft denotes any combination 
of Clifford algebra elements and color group generators with a combined index J which is contracted. 
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C. Spurion analysis 



So, the leading term in the Symanzik action, 6 is just the continuum action of partially 
quenched QCD. In the m — > limit it is invariant under the flavor symmetry group Gpq in 
eq. (|23|) . which is larger than Gm [eq. (|24j)]. the symmetry group of the underlying lattice 
action. For a sufficiently small a (and m), So determines the spontaneous symmetry breaking 
pattern and the symmetry properties of the Nambu-Goldstone particles in the theory It 
follows that the mixed theory contains the same set of light particles as partially quenched 
QCD. 

For the construction of the chiral effective theory we introduce spurion fields that make 
the entire Symanzik action invariant under Gpq. Notice that all the operators proportional 
to a and a 2 break Gpq, the flavor symmetry of the leading term. This is obvious for operators 
which appear with sea quark fields only, like the dimension 5 operators. However, even if 
an operator appears "symmetrically" in if)$ an d ipv, as i n eq. (EBj) ; it still breaks Gpq. To 
illustrate this point let us consider any of the bilinear terms, suppressing all 7 matrices and 
color group generators. Any bilinear that is invariant under all rotations of Gpq must have 
the flavor structure = ipgips + VvVv- in general, though, ip s ips and ipyipv will not 
appear in the Symanzik action with equal coefficients, and therefore will not be invariant 
under transformations in Gpq that mix the sea and valence sectors. 

As before we begin the construction of the chiral Lagrangian by listing the representative 
spurions required at each order in a to make the Symanzik action invariant. Shown are 
the transformation properties of the spurions under chiral transformations in Gpq and the 
constant structures which the spurion fields are assigned to in the end. Since different 
operators appear in the sea and valence sector it is convenient to introduce the projection 
operators 

P s = diag(/ 5) 0) , P v = diag(0, I v ) , (30) 

where Is denotes the Nf x Nf identity matrix in the sea sector, and Iy the 2Ny x 2Ny 
identity matrix in the space of valence quarks and ghosts (recall that tpy includes both 

6 As before, we assume that the sea quark mass renormalization effectively accounts for S-\- see eq. <|11[) . 
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valence quarks and ghosts). 

0(a°) : M -> LMR\ M ] -> RM^lJ, (31) 

M = Mq — m — diag(m5, m y ). 

0(a) : A -> LAR\ A ] -> P^L 1 , (32) 

A = = aP s . 

The last spurion arises from the sea sector symmetry breaking terms at 0(a). 

All the quark bilinears at 0(a 2 ) couple fields with the same chirality. Since there are 
bilinears for both sea and valence fields we obtain the following spurions: 

0(a 2 ), bilinears : B -> LBL\ C -> RCR\ (33) 

B , C E {a 2 P s ,a 2 P v }. 

We can distinguish two types of 4-quark operators. The first type is made of bilinears 
which couple fields of the same chirality only. These operators appear with only sea or 
valence fields as well as in the "mixed" form 0(ips, Vv)- The remaining 4-quark operators, 
which couple fields with opposite chirality, appear only with sea quarks. We therefore 
introduce the following spurions: 

0(a 2 ), 4-quark operators : 

D = Di <g> D 2 -> LD\l) ® LD 2 L\ E = E l ®E 2 ^ RE^ ® RE 2 R\ 
F = F 1 ®F 2 ^ LF X I) ® RF 2 R\ G = d <g> G 2 -f PG^P 1 ® LG 2 L\ 
Do, E , Fq, Go E {a 2 P s ® P 5 , a 2 P 5 ® Py, « 2 ^y ® (34) 

H = Hi® H 2 ^ LH^ <g> LiJsPt, P f = ffj ® #| -> PP 1 t L t ® PPf^, 
J = P. ® j| -> LJi^ ® Rj\L\ J f = 4 ® J 2 -> i?4^ ® J^P 1 , 
#o = PfJ = Jo = 4 = a 2 P s ® P?- (35) 
Squaring the spurions of 0(a) does not lead to any new spurions. 



D. Chiral Lagrangian 



The chiral Lag rang ian for the mixed action theory including the cut-off effects linear in a 
is derived in Ref. 16|. Terms of 0(a 2 ) are constructed from the spurions in eqs. (j33H35|) . It 
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is easily checked that the spurions B, C, D and E lead necessarily to operators higher than 
0(a 2 ) [at least 0(p 2 a 2 , ma 2 )], so we can ignore them. From the other spurions we obtain the 
following independent operators (and their hermitian conjugates): 

(F 1 EF 2 ^) - a 2 <r 3 Sr 3 St) , (36) 

(HjfHttf) - a 2 (P s £tp s £t) ? (37) 

(Hrf) (H£$) - a 2 (P s Et) ( PflE t) ) (38) 

< JJ?) ( 4s\ - a 2 <P 5 St) (P S E) . (39) 



For eq. we use the fact that P5 = |(/ + r 3 ) and Py = §(/ — t 3 ), with r 3 = diag(J s , — Iy). 
When assigning F\ t 2 = (I ± r 3 ) and expanding, the fields £ and E^ are next to each other 
and cancel whenever the identity matrix is inserted, so the only non-trivial operator is the 
one shown in eq. ([3*B|l. 

We conclude that for the mixed action theory with Wilson sea and Ginsparg- Wilson 
valence quarks the terms of 0(a 2 ) in the chiral Lagrangian are 

£ [a 2 ] = -a 2 W M <r 3 Sr 3 St) (40) 
- a 2 Ws(P s ^ + SP S ) 2 - a 2 W;(P s ^ - SP 5 > 2 - a 2 W^(P s ^P s ^ + EP s EP<j}. 

The parameters rh and a are defined as in the unquenched case in eq. ([20)1 . Note that the 
projector P$ in the last three terms implies that these operators involve only the sea-sea 
block of E. 

The final result, including the terms from Ref. [lfj], reads 

L x = £ (d^Zd^) - ^ (m£t + Em) - c£ {P s tf + P 5 S) 

- L, {d^d^f - L 2 (^E^E+) (0 M £&£t) - L 3 ((^E^Et) 2 ) 
+ L 4 (<9 /1 E(9 M £ t ) (mE f + Em) + al^ 4 (9 M E9 M E t ) (PsS 1 " + SP 5 ) 
+ L 5 (^E9 M E+(mEt + Em)) + alf 5 (d^E^P^ + SP 5 )) 

- L 6 (m£ f + Sm) 2 - aW 6 <mS f + Em) <P 5 S f + EP 5 ) 

- L 7 (mS f - Sm) 2 - aVT 7 (m£ f - Em) (P 5 E f - EP 5 ) 

EmEm) - aW 8 (mE 1 P s E f + EPgEra) 
+£[a 2 ] + higher order terms. (41) 
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The chiral Lagrangian for the mixed action theory at 0(a 2 ) has 4 terms while there are 
only 3 terms at this order in the chiral Lagrangian for the Wilson action. The reason that 
the mixed theory has an additional operator (and consequently an additional unknown low- 
energy constant multiplying it) is its reduced symmetry group, Gm in eq. (|24|). compared 
to Gpq in eq. (|23|) . The use of different Dirac operators for sea and valence quarks for- 
bids transformations between the sea and valence sectors, and allows the additional term 
(r 3 Sr 3 St) in eq. flUD. 

The presence of more terms in the Lagrangian does not entail that chiral expressions for 
all observables in the mixed theory depend on more free parameters than in ;^PT for the 
Wilson action. By definition, the correlation functions measured in numerical simulations 
involve operators that are made of valence quarks only, and the enhanced chiral symmetry of 
the Ginsparg- Wilson fields plays an important role in that sector. The chiral symmetry leads 
to constraints on operators in the Symanzik action that contain valence fields, and ultimately 
it restricts and simplifies the form of chiral expressions for valence quark observables. This 
can already be seen by considering the terms in L [a 2 ] with coefficients Wq, and Wg [see 
eq. (|4()j) ]. These terms depend only on the sea-sea block of S. This entails that all the 
multi-pion interaction vertices obtained from these terms necessarily contain some mesons 
with at least a single sea quark in them. Consequently, these terms cannot contribute at 
tree level to any expectation value of operators made entirely out of valence fields. This is 
easily understood: the W terms arise from the breaking of chiral symmetry in the sea sector 
by the Wilson term, and this breaking is communicated to the valence sector through loop 
effects only. A more concrete demonstration of this point is provided by the calculation of 
the pseudo scalar valence-valence meson mass in the next section. 

IV. APPLICATION 

We conclude our analysis of the chiral effective theories for the Wilson action and the 
mixed action theory with an explicit calculation of the light meson masses. Before presenting 
the calculations, however, a discussion of the chiral power counting is appropriate. 
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A. Power counting 



%PT reproduces low-momentum correlation functions of the underlying theory, provided 
that the typical momentum p and the mass of the Nambu-Goldstone boson M^gb are 
sufficiently small, p < A x and M NGB <C A x . The standard convention is to consider p and 
Mngb as formally of the same order, and take a single expansion parameter e ~ M]^ GB / A x ~ 
p 2 /A^. Thus, a typical next-to-leading order (1-loop) expression for a correlation function 
in %PT has the structure 



In some cases of interest the momentum scale and the Nambu-Goldstone boson mass 
are significantly different, p <C Mmgb for instance. In such a case one could treat the two 
dimensionless parameters separately, and introduce another expansion parameter, p/M NGB . 
However, so long as both Mj fGB /A^, and p 2 7 'A 2 are sufficiently small, eq. still holds. 
Consequently, a reasonable approach in the case that p and M^gb are very different is to 
take eq. (|42J1 and to ignore (or not calculate) terms that are smaller than the error associated 
with the larger expansion parameter. 

In the case of for lattice theories there are two possible sources of explicit chiral 
symmetry breaking: the quark masses and the lattice spacing. Consequently, the mass of 
the pseudo Nambu-Goldstone boson is given by M^ GB /A 2 ~ m/A x + aA x . The discussion 
of the previous paragraph applies here as well: we can take e ~ 7 'A 2 , ~ m/A x ~ aA x 
and eq. (|4"2*j) (properly extended) still holds. So long as the largest of these parameters is 
sufficiently small this is a consistent power-counting scheme, and Eq. (}4"2*j) is applicable even 
when some of the dimensionless parameters are significantly smaller than the others. This 
is the power-counting which is used in organizing the terms in the Lagrangians in eq. (|22|) 
and eq. (JUJ). 

A different power counting scheme does need to be employed in some cases. To illustrate 
this we consider a realistic example: for some fermion actions there are no discretization 
effects at 0(a). This is the case, for example, for non-perturbatively 0(a) improved Wilson 



C 



Clo + Cnlo + • • • j 



(42) 
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fermions. If, in addition, the lattice spacing in a simulation is large such that a 2 A 2 ~ m/A x , 
an expansion in two parameters may be required, and the leading order contributions are 
0(p 2 /A 2 ,m/A x ,a 2 A 2 ). 



B. Pseudo scalar meson masses 

We now turn to the calculation of the pseudo scalar meson masses. As in Refs. Ml and 
, we only consider mesons with different valence flavor indices (A ^ B). For the partially 
quenched Wilson action we find 

m ab = (™vai + a) + 9 (m Va i + a) [m Va i - m ScSu + (2m Va i - m Sca + a) In (m Va i + a)] 

g 

--^(m Va i + a) [N f (L 4 m Sca + W A a) + L 5 m V ai + W 5 a) 

+ ^-[2L 6 m V ai^Sea + W 6 (fh Yal + m SeSL )a + 2W^a 2 ] 
16 

+ — [L 8 m 2 al + W 8 m V aia + W$ 2 } + 0(e 3 ), (43) 

where Nf is the number of sea quark flavors. 

Next, we consider the mixed action theory. As explained in section IIII1 the terms in 
the Lagrangian at 0(a 2 ) with coefficients W§, W!j and W$ [see eq. pOjl] cannot contribute 
at tree level to any expectation value of operators made entirely out of valence fields. In 
particular, at next-to-leading order they do not enter the self energy of mesons made out 
of valence quarks, and thus do not contribute to their mass. In addition, the leading order 
meson mass (M| b )lo i n the mixed action theory is independent of a, which means that 
contributions at 0(e 2 ) of the type «(A^vgb)lo (from the terms proportional to W4 and W 5 
) are only linear in a. Thus, the only possible source for a 2 corrections to the mass comes 
from the term in SL, [a 2 ] with coefficient Wm- However, an explicit calculation shows that 
this contribution vanishes and there are no 0(a 2 ) corrections to the pseudo scalar meson 

n 

mass. The expression for Mab is therefore the same as in Ref. [16J, which we quote here for 
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completeness: 7 

M AB = "^Val + T^TT-FTT^Val [m Va J ~ m Soa ~ a + (2m V al - "^Soa - a) In (m V ai)] (44) 

-^m Val [(i5 - 2L 8 )m Va i + N f (L 4 - 2L 6 )m S ea + iV/(W 4 - W & )a] + 0(e 3 ). 

The fact that there are no 0(a 2 ) contributions at next-to-leading order is not as surprising 
as one might think at first. Only the valence quark mass term breaks the chiral symmetry 
for Ginsparg- Wilson fermions. Hence the pseudo scalar meson mass is proportional to the 
quark mass and vanishes in the limit my a i — > 0. It follows that any lattice contribution to 
M\ B is suppressed by at least one factor of m w , and the largest lattice correction quadratic 
in the lattice spacing is of 0(my a /a 2 ). Note that this higher order term becomes the leading 
discretization effect in the meson mass if an 0(a) improved Wilson action is used for the 
sea quarks. This example illustrates the beneficial properties of Ginsparg- Wilson fermions 
which are preserved even in the presence of a "non- Ginsparg- Wilson" sea sector. 



V. SUMMARY 



In the previous sections we presented the chiral Lagrangian for two lattice theories: one 
with Wilson fermions and the other with Wilson sea fermions and Ginsparg- Wilson valence 
fermions. One consequence of the analysis is that corrections to the low-energy constants 
of continuum (coming from symmetry-conserving discretization effects) are of 0(a 2 ). 
Since the coefficients in the chiral Lagrangian themselves multiply terms of 0(p 2 ) (Bq and 
/) and 0(p 4 ) (Gasser-Leutwyler coefficients), such effects can only be detected by measuring 
observables at the accuracy of 0(a 2 p 2 ) and 0(a 2 p 4 ), respectively. Another important dis- 
cretization effect that enters the Symanzik action at 0(a 2 ) is the breaking of 0(4) rotational 
invariance. An 0(4) breaking term in the chiral Lagrangian, however, must contain at least 
four derivatives, so it is a higher order term as well (at least 0{a 2 p i )). 

The main purpose of constructing chiral effective theories for lattice actions is to capture 
discretization effects analytically and to guide the chiral extrapolations of numerical lattice 
data. This is achieved by the explicit a dependence of observables that can be calculated in 
these effective theories. In particular, the chiral Lagrangian is sufficient for the determina- 

7 In Ref. the number of flavors Nf was set to three. 
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tion of the pseudo scalar meson masses. For the calculation of matrix elements like f n an 
additional a dependence coming from the effective continuum operators needs to be taken 
into account, but no conceptual difficulties are expected to arise in this step. 

There is a more subtle cut-off dependence that is not explicit in the Symanzik action. 
All the unknown coefficients in the Symanzik action implicitly include short- distance effects 
that make them a dependent. An important example is csw{a) for which the a dependence 
has been calculated non-perturbatively in quenched and unquenched simulations. For the 
chiral Lagrangian this results in an implicit a dependence of the low-energy constants 
The existence of a well-defined continuum limit implies that all the parameters of continuum 
xPT, such as the Gasser-Leutwyler coefficients, have a leading, a-independent part. The 
other coefficients in the Lagrangian, loosely referred to as the "W's" , are expected to show 
a weak, logarithmic a dependence. 

From a practical point of view there are several ways to approach this issue. One option 
is not to vary a. For a given lattice spacing a one fits the chiral forms by only varying the 
quark masses. Note that even if a is not varied, the inclusion of the discretization effects 
in the chiral expressions, particularly in the chiral logarithms, is more accurate than simply 
using the continuum expressions. From the fits one extracts values for the coefficients in the 
Lagrangian, including the W's. Applying this procedure again, independently and for lattice 
data with different lattice spacings, these parameters are allowed to vary with a. It should 
be verified that the values obtained in this way for the continuum low-energy constants do 
not exhibit an a dependence beyond the error expected at the order of the calculation. It 
might be the case that the a dependence of the W's is so slow that they do not change 
much over the range of lattice spacings simulated. In that case a simultaneous fit in a 
and m might be appropriate. More generally, a simultaneous fit can be used when the a 
dependence of the W's is known. In particular, it is consistent to treat the W's that enter the 
chiral Lagrangian at 0(a) as being proportional to c$w- This step simply amounts to using 
equations of motion in the Symanzik action. This way, since csw( a ) is numerically known, 
one might have control over all the a dependence in the chiral Lagrangian at 0(a), pushing 
the unknown a dependence to 0(a 2 ). This is "automatically" done in O(a)-improved lattice 
simulations. 

All the qualifications of the previous paragraphs notwithstanding, chiral perturbation 
theory for lattice actions provides a better understanding of the relation between lattice 
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observables and their continuum counterparts. It is encouraging that at 0(a 2 ) only a few 
new low-energy constants are needed. Thus x?T is still predictive at this order and it is 
likely to play an important role in the extraction of quantitative predictions of QCD from 
numerical simulations. 
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APPENDIX A: FLAVOR, COLOR AND DIRAC STRUCTURE OF 4-QUARK 
OPERATORS IN THE SYMANZIK ACTION 

In this section we discuss 4-quark operators in the Symanzik action that are invariant 
under the vector flavor symmetry group SU(Nf)y, the color gauge group SU(N C ), the 
hypercubic transformations, parity, and charge conjugation. 

It is convenient to label the quark fields ^\ ■ Considering first the flavor 

group, we write the most general term (summation over repeated indices is assumed) 

There are only two possibilities for C (up to a multiplicative constant) which make this term 
invariant: 

Cili2i3«4 ^1112^1314) and Ci 1 i 2 i a i 4 ^i^i^i^- (A2) 

These correspond (up to a sign from the interchange of the Grassman fields) to 

^VV/Vf and ^V^fVf. (A3) 

At this point we are free to redefine the labels on the quark fields in the second term by 
exchanging the second and fourth indices. This way we only need to consider the first 
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invariant in the last equation. From this point on the order of the fields will remain fixed, so 
the labels of the fields can be dropped, and the trivial flavor contractions will be suppressed. 

The same analysis holds for the color structure. The difference is that now we have 
already exhausted the freedom to reshuffle and relabel the fields - they are distinguishable 
by their flavor indices - and so there are two genuinely different invariant operators: 

^a^a^ b ^b, and ip a i/; b ip b ip a . (A4) 

We find it convenient to "untwist" the color indices in the second term using the Fierz rule 

SaJbd = ^abtcd + KbKd, (A5) 

where the t e are the generators of the color group in the fundamental representation. The 
possible terms can now be written as 

ifyijrif), and ijt a "ipi)t a "ip, (A6) 

where the contraction of color indices is straightforward and can also be suppressed. 8 

Finally, we take into account the Dirac structure. To maintain the hypercubic symmetry 
and parity invariance, the space-time indices must be contracted in pairs and 75 matrices 
must appear in pairs. One set of invariant terms can be obtained by adding a Dirac structure 
to the terms in eq. (jA6|) in the following way: 

^T A ijj^T A ilj and ^V A t a ^Y A t a ^, ( A7 ) 

where ijjY A ij) can be a scalar, pseudoscalar, vector, pseudovector, or a tensor, with A denoting 
the appropriate space-time indices. In addition, as in the cases of color and flavor, it is also 
possible to have the Dirac matrices connect the first and fourth fields, and the second and 
the third. These operators, however, are linearly dependent on the previous terms, because 
of the identity 

V A V B _ \^ T^AB r C pfl 
1 aS L 7/3 — CD a/3 j5i 

C,D 

K A % = — Ti[T A T D T B T c ). (A8) 
16 



In fact, the color structure is completely inconsequential in the construction of the chiral Lagrangian. 
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This identity holds for any pair of Clifford algebra elements and not only for the case A = B 
in which we are interested. 

This completes the derivation of the list of 4-quark operators in eq. (jl(J|) . There are several 
equivalent sets of operators. A different path leads to the list of operators that appear in 
Ref . [3] : starting with the color structure, one considers the invariants of eq. (jA3|) but with 
color indices instead of the flavor ones. Fierz rules can be used to replace the identity 
matrices with color generators that are either "straight" (connecting the first and second 
fields, and the third and fourth) or "crossed". As was done above with respect to flavor, 
it is possible to choose a convention in which all the color generators are straight (reorder 
the fields). Thus the only invariant is ijjt a tjjtjjt a tp. Once that convention is fixed, one again 
faces the possibility of crossed Dirac and flavor indices. The Dirac matrices are straightened 
in the same manner as above. Finally, using Fierz rules for the flavor group, one can also 
eliminate the crossed Kronecker deltas at the price of introducing terms with flavor group 
generators (3\ The final set of invariants is ipT A t a il)il)T A t a ii and ifjT A t a j3 i ipiljT A t a l3 i i/j. 

APPENDIX B: REDUNDANCY OF SPURIONS 

We note the following fact: If A and B are two spurions which are of the same order 
in the (m, a) power counting, transform in the same manner, and have a similar "original" 
structure, Bq = kAo where k carries no indices, then one can use only one of them to 
construct the chiral action. The reason is the following. If f(A) is an operator in the chiral 
Lagrangian which contains A, then f(B) is also an allowed term - because of the assumption 
that both spurions transform in the same way. Since the spurions transform linearly, the 
relation between the constants A and B Q also holds for the spurions. Assuming a power 
expansion in the spurions, this leads to f{B) = k n f(A). Recalling that each operator in the 
chiral Lagrangian appears with an unknown coefficient, we have in the Lagrangian 

K x f{A) + K 2 f(B) = (Kt + K 2 k n )f(A). (Bl) 

Since neither K\ nor K 2 are known (and in most cases neither is k), this is equivalent to 
considering only a single term in the Lagrangian, Kf(A), which we would have written 
anyway if we had considered only the first spurion. 
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Example: At 0(a) the Symanzik Lagrangian contains the terms 

acii) L D^D^ R + ac 2 ^ L i(T flu F^ R . (B2) 

To make these terms invariant one can introduce two spurions A and B which are flavor 
matrices that transform as A — > LAR\ B — > LBR) . Both are 0(a), and their constant 
values are A = ac\I, B Q = CLC2I (here / is the flavor identity matrix). With these 
spurions it is possible to construct the following invariant terms in the chiral Lagrangian (at 
0(a)): 

K 1 {AY I ] ) + K 2 {BT?) (B3) 
but after setting the spurions to their constant values we obtain only a single term 

K x aci (£ f ) + K 2 ac 2 <£ f ) = aK' <£ f ) (B4) 
which we would have writing down even if we had kept only A in the analysis. 
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